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ABSTRACT 

Sun-like stellar oscillations are excited by turbulent convection and have been discovered in some 
500 main sequence and sub-giant stars and in more than 12,000 red giant stars. When such stars 
are near gravitational wave sources, low-order quadrupole acoustic modes are also excited above 
the experimental threshold of detectability, and they can be observed, in principle, in the acoustic 
spectra of these stars. Such stars form a set of natural detectors to search for gravitational waves 
over a large spectral frequency range, from 10~^ Hz to 10“^ Hz. In particular, these stars can 
probe the 10“^ Hz - 10“^ Hz spectral window which cannot be probed by current conventional 
gravitational wave detectors, such as SKA and eLISA. The PLATO stellar seismic mission will achieve 
photospheric velocity amplitude accuracy of cm/s. Eor a gravitational wave search, we will need to 
achieve accuracies of the order of 10“^cm/s, i.e., at least one generation beyond PLATO. However, 
we have found that multi-body stellar systems have the ideal setup for this type of gravitational wave 
search. This is the case for triple stellar systems formed by a compact binary and an oscillating star. 
Continuous monitoring of the oscillation spectra of these stars to a distance of up to a kpc could lead 
to the discovery of gravitational waves originating in our galaxy or even elsewhere in the universe. 
Moreover, unlike experimental detectors, this observational network of stars will allow us to study the 
progression of gravitational waves throughout space. 

Subject headings: asteroseismology - binaries: close - binaries: general - gravitational waves - stars: 
general - Sun: helioseismology 

1. INTRODUCTION 

Astronomers had not expected that stars, other than 
the Sun, could provide us with the observed wealth of 
high quality seismic data, surpassing in diversity and 
quantity the data gathered from the Sun itself. Previ¬ 
ously, the COLE experiment aboard the Solar and He¬ 
liospheric Observatory mission Gabriel et al. (1995) mea¬ 
sured the global acoustic oscillations of the Sun from 
a spatially unresolved entire solar disk of Doppler ve¬ 
locity data (Turck-Chieze et al. 2004). Presently, fol¬ 
lowing an identical strategy, the COROT and Kepler 
space missions (Baglin et al. 2006; Gilliland et al. 2010), 
by measuring the light integrated from the entire visi¬ 
ble star’s surface, have searched for global oscillations 
in more than 150,000 main sequence, subgiant and red 
giant stars (Borucki et al. 2009; Verner et al. 2011). A 
sample of more than five hundred low mass subgiants and 
main sequence stars were discovered to have rich Sun-like 
acoustic oscillation spectra filled with tens of velocity am¬ 
plitude peaks (Chaplin et al. 2014). Similar oscillations 
have been found in more than 12,500 K-G giant stars. 

Although their acoustic spectra are very distinct from 
those of main sequence stars, as a consequence of their 
quite distinct structure (Hekker et al. 2011), these oscil¬ 
lations still qualify as Sun-like oscillations. As in the Sun, 
these oscillations are excited by turbulent convection and 
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intrinsically damped by the radiation and the turbulence 
of the outer layers of the star. Such seismic observational 
surveys are able to follow in a systematic and continu¬ 
ous manner the pulsation spectra of the Sun and stars 
in the solar neighbourhood within a range of up to one 
thousand parsecs distance. These surveys observed fields 
of stars in the galactic plane, as well as in a few direc¬ 
tions above and below it (Miglio et al. 2012), as in the 
Kepler mission which has been observing stars 13° above 
the galactic plane in the Cygnus region. This local grav¬ 
itational system of celestial bodies made of many stars 
of different sizes and masses, constituted of single stars, 
binaries or multi-stellar systems, as well as a few stellar 
clusters (Corsaro et al. 2012), is a natural network of de¬ 
tectors for gravitational radiation. In this work we argue 
that this natural network of stars will allow us to search 
for gravitational wave imprints in the oscillation spectra 
of single stars, or for contemporaneous signatures of the 
same gravitational wave event on the spectra of two or 
more close stars. It will also allow us to follow the pro¬ 
gression of a gravitational wave through space, across the 
spectra of several stars - when it passes through a stellar 
cluster - and to monitor its progression as it approaches 
the Earth. 

The most likely source of gravitational waves able to 
stimulate quadrupole modes in Sun-like stars are short- 
period binaries of two compact objects such as white 
dwarfs, or even neutron stars and massive black holes. 
Low mass binaries are the leading gravitational source 
candidates, despite having relatively weak gravitational 
wave emission, as they are numerous and are located in 
close proximity to the Solar System (Amaro-Seoane P. 
et al. 2013). More than 50 ultra-compact binaries, with 
a size of about a fraction of the solar radius and a period 
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shorter than one hour, have been discovered at a distance 
between 50 and 700 parsec, an example being the cata¬ 
clysmic variable star AM CVn located at a distance of 
606 parsecs from Earth (Roelofs et al. 2007). AM CVn 
is a binary system where a white dwarf accretes matter 
from a companion star, leading to the formation of an ac¬ 
cretion disc with continuous strong emission in UV and 
X-rays, and with occasional outbursts. If seismic surveys 
are set to observe and monitor stars near these binaries, 
this could lead to the discovery of gravitational waves. 

The idea that gravitational radiation could excite the 
normal modes of vibration of celestial bodies such as 
the Earth and the Sun was originally discussed by 
Dyson (Dyson 1969) and many other papers have fol¬ 
lowed up this idea. Most recently McKernan et al. (2014) 
have estimated the gravitational radiation that is ab¬ 
sorbed by stars and Siegel & Roth (2011) are among 
others to suggest the impact of gravitational waves on 
solar oscillations. Gravitational wave detection through 
stars and resonant mass detectors (Aguiar et al. 2006; 
Gottardi 2007) both work based in a similar principle. 
In the latter case, the detection is done by accurately 
measuring the tiny variation of the detector size when 
one or more modes are excited by a passing gravitational 
wave (Sathyaprakash & Schutz 2009). In stars detection 
is feasible by monitoring the variations of velocity of the 
modes at the surface. Some aspects of the analysis for 
stars are similar to the case of spherical resonant mass 
detectors. 

2. THE ACOUSTIC SPECTRA OF THE SUN AND STARS 

Stars, like musical instruments, vibrate in a multitude 
of eigenmodes. The discrete sequence of frequencies for 
stellar oscillations can be labelled by two independent 
integers: the degree of the mode I, which is the degree of 
the spherical harmonic related with the horizontal eigen¬ 
function; and the order of the mode n which measures 
the number of nodes of the radial eigenfunction along 
the radius. Eor each value I, there is a sequence of reso¬ 
nant acoustic modes that are labelled with n. The latter 
ones correspond to the fundamental tone and overtones 
of a musical instrument. Whole-disk observations of stars 
through the seismic space missions detect only the large 
luminosity variations in the stellar surface; therefore, 
these observations are only sensitive to the lowest values 
of I (I ^ 3). The angular frequency of such low-degree 
acoustic modes satisfies 

CUlM = (r + 1/2 I -h (Xo)d) + Cn (1) 

where N is a subscript that defines a specific eigenmode 
N = (n, 1), ao and d) are constants and Cn is a second 
order term that can be neglected when n is large (Lopes 
2001). The constant d) is relates do Vq by Vq = Cd/2n. 
This last quantity is known as the large separation. The 
frequencies with the same degree are separated apart by 

Vq. This quantity is equal to 1/ ^2 dr/Vg (r)^, the 

time taken for the sound wave to travel with a speed 
Vs(r) from the surface to the center of the star and re¬ 
turn. According to equation (1) for cum, Vq is equal to 
(H)(n,i) ~ ]^))/27r. The previous equation with mi¬ 

nor adjustments has been shown to be valid for many 
main sequence and red giant stars (Mosser et al. 2011; 
Gorsaro et al. 2012). Table 1 lists the frequencies of 


the quadrupole mode as measured by the GOLE ex¬ 
periment. When no observational data is available, we 
show in italics the predicted values for the standard so¬ 
lar model (Lopes & Turck-Ghieze 2013). One can notice 
that for frequencies of the acoustic modes of low degree, 
the disagreement between theory and observation is at 
most of the order of a few percent (see Turck-Ghieze 
& Lopes 2012, and references therein). Table 1 lists the 
frequencies for the quadrupole modes. 

Recent observations have shown that many main- 
sequence, subgiant and red giant stars have a spectrum of 
acoustic oscillations identical to the Sun, with only minor 
differences (Mosser et al. 2011; Chaplin et al. 2014). The 
properties of a star’s spectrum are characterized by three 
main quantities: the large separation, the frequency at 
which the amplitude of the spectrum of oscillations is 
maximum, and the effective temperature of the star. At 
present, when these quantities are available, they provide 
the most reliable method for determining the mass and 
radius of a star (Chaplin & Miglio 2013). It has been 
shown from observational data, that homology ratios for 
mass, radius and effective temperature hold between the 
Sun and these stars. In particular, the large separation 
scales as: 



where M. and R are the mass and radius of the star. 
Equally, cDQ/27r, Mq and R© are the solar equivalent 
quantities. In the case of the Sun, the mean large sep¬ 
aration (x)/2n is of the order of 135 p Hz (as computed 
from table 1). 

3. INTERACTION BETWEEN GRAVITATIONAL WAVES 
AND ACOUSTIC MODES 

Within the framework of general relativity, far from the 
source of gravitational radiation, the space-time metric 
tensor is distorted relative to flat spacetime (Minkowski) 
value by a very small spatial component, hij. In a 
Galilean coordinate frame whose origin coincides with 
the center of the star, the stellar material experiences a 
force proportional to hij. Thus, the quadrupole modes 
of vibration of the star will be excited by gravitational 
waves when the frequency of the incoming waves is close 
to the eigenfrequency of the modes. In this case, the 
time variation of the amplitude of the mode of vibration 
is described by a harmonic damped oscillator with an ex¬ 
citation source proportional to hij. As usual, we express 
the tensor hij as the sum of spherical components hra, 
for which the m (azimuthal order) is an integer such that 
|m| ^ 1. The dynamics of general relativity implies that 
only non-radial modes of degree larger that two can be 
excited. Of these, the forcing of the quadrupole modes is 
normally the greatest. Actually, the differential rotation 
in the Sun and Sun-like stars can also split the cum (equa¬ 
tion 1) leading to a subset of 2m + 1 frequencies for each 
I, which in the case of quadrupole modes corresponds 
to five values differing between them only by a few mi- 
croHz. Nevertheless, as we are only concerned about the 
amplitude of the modes, in this study the solution will be 
found for a generic value of m. As we restrict our atten¬ 
tion to quadrupole acoustic modes of order n with fixed 
(undistinguished) m for which |m| ^ 2, in the remainder 
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TABLE 1 

Quadrupole acoustic modes 
Observational data and Standard solar model 


n 

Freq. ^ 

n b 

On 

Xn 

Ln 

Mn 

In 

Qn 

Vn " 

"T* d 

'n 

V ® 

V n,u 


(M-Hz) 

(M-Hz) 

(no-dim) 

(cm) 

(Mo) 

(cm“^ Hz) 

(no-dim) 

(cm ) 

(no-dim) 

(cm ) 




xl 0 -"‘ 

xlO’’ 

xlO-^ 


Xl0+^ 

XH _20 

xlO-^ 

XH _20 

f 

347.10 

2.9 10-'' 

-6.7432 

2.347 

0.5854 

94.931 

38 

1.410-8 

2.7 

7.1 10-8 

Pi 

382.26 

2.9 10-'' 

-11.038 

3.841 

1.1148 

98.315 

4.1 

2.610-1 

23 

4.0 10-8 

P2 

514.48 

1.6 10-" 

±2.1193 

0.737 

0.1594 

0.9671 

1.0 

1.710-8 

72 

4.6 10-° 

P3 

664.06 

7.4 10 -" 

-0.6286 

0.219 

0.0466 

0.0424 

0.28 

1.810-° 

210 

8.2 10 - 1 ° 

P4 

811.33 

2.6 10-^ 

+0.2133 

0.074 

0.0154 

0.0025 

0.10 

2.610-1° 

506 

1.9 10-1° 




xl 0 “^ 

xlO^ 

xlO-^ 


xlO^ 


xlO^ 


P5 

959.23 

7.9 10-^ 

-8.2377 

2.867 

1.0733 

2 10 ^"^ 

3.8 

4.710-11 

1.1 

(5.0 10-11) 

P6 

1104.28 

2.1 10-" 

±3.4804 

1.211 

0.4537 

2 lO-'^ 

1.7 

9.810-12 

2.3 

5.5 10 - 11 ) 

P7 

1249.78 

5.3 10-" 

-1.5051 

0.524 

0.1968 

2 10 “® 

0.74 

2 . 110-12 

4.7 

5.6 10-12) 

P8 

1394.68 ± 0.01 

0.01 

±0.6836 

0.238 

0.0880 

2 10 “^ 

0.33 

4.910-1® 

9.4 

5.5 10 - 12 ) 

P9 

1535.865 ± 0.006 

0.04 

-0.3109 

0.108 

0.0392 

3 10-® 

0.12 

8.910-11 

24 

(4.4 10-1®) 




xlO-^ 

xlO^ 

xlO-^ 


xlO^ 


xlO^ 


PlO 

1674.534 ± 0.013 

0.08 

+ 14.946 

1.082 

18.647 

3 10-® 

6.7 

2.610-11 

4.1 

— 

Pll 

1810.349 ±0.015 

0.10 

-7.8242 

0.520 

10.06 

6 10 - 1 ° 

4.1 

9.010-18 

6.2 

— 

P 12 

1945.800 ± 0.02 

0.14 

+4.3862 

0.272 

6.106 

2 10 - 1 ° 

2.9 

3.810-15 

8.3 

— 

Pl3 

2082.150 ±0.02 

0.21 

-2.5981 

0.153 

4.0176 

3 10-11 

2.1 

1.710-15 

10.9 

— 

Pl4 

2217.69 ± 0.03 

0.31 

+ 1.5564 

0.054 

2.7101 

9 10-12 

1.9 

1.010-15 

11.4 

— 

Pl5 

2352.29 ± 0.03 

0.40 

-0.9562 

0.033 

1.9063 

3 10-12 

1.8 

6.410-18 

11.1 

— 

PI 6 

2485.86 ± 0.03 

0.44 

+0.6204 

0.022 

1.4536 

1 10-12 

1.7 

4.210-18 

11.1 

— 

Pit 

2619.64 ± 0.04 

0.53 

-0.4180 

0.014 

1.1950 

4 10-1® 

1.5 

2.610-18 

12.2 

— 

P 18 

2754.39 ± 0.04 

0.57 

+0.2908 

0.010 

1.0383 

2 10 - 1 ® 

1.5 

1.810-18 

12.0 

— 


^The observational frequency corresponds to (X)^/2n. The table is obtained from a compilation made by Turck-Chieze & Lopes ( 2012 ), 
after the observations of Bertello et al. (2000); Garcia et al. (2001); Turck-Chieze et al. (2004); Jimenez & Garcia (2009). The frequencies 
in italic correspond to theoretical predictions for the current standard solar model as in reference Lopes & Turck-Chieze (2013). 

^The damping rates are interpolated from a rj observational table of averaged values obtained for all global modes with 1^3 (Chaplin 
et al. 1997). These observational results are consistent with the rj values of dipole modes measured by Baudin et al. (2005). The theoretical 
values of damping rates are from Houdek et al. (1999); Grigahcene et al. (2005); Belkacem et al. (2013). 

^The photospheric velocity is computed for a strain of H 0 = H_ 2 o 10“^° with 1l_2o ~ 1- 

‘^This value is computed for the compact binary system AM CVn (see Figure 1 for details). 

®This photospheric velocity corresponds to the unsaturated limit T* ^ 1 (or Tf ^ Tn). The values within (•) are indicated for reference 
only. 


of this article the fiducial mode will be represented as N 
or simply by n, meaning N = (n, 2, ”m”). 

The strength by which the quadrupole modes of the 
star are stimulated by gravitational radiation depends on 
the absorption cross-section for gravitational radiation 
or its integrated value (in frequency), the quality factor, 
and the amplitude of the root mean square velocity at 
the surface of the star, also known as the photospheric 
velocity. The theoretical calculation of these quantities is 
computed in an identical manner to gravitational wave 
detectors. Accordingly, the impact of a plus polarized 
monochromatic gravitational wave, such as h* cos (cut), 
with a strain and frequency cu, averaged over several 
cycles is estimated as follows: 

- First, the absorption cross-section for gravitational 
waves by the star, crabs(cu), is defined by expressing the 
balance between the amount of energy Eabs which is ab¬ 
sorbed by the star and the amount of incident energy Ein 
on the star’s surface: 


dEabs 

dt 


0‘abs(cu) 


dEia 

dAdt' 


(3) 


where dEin/dAdt is the energy arriving per unit of time, 
per unit area. 

- Second, the average dEin/dAdt value for the case 
of a monochromatic gravitational wave is equal to 
c^H^cu^/(327rG) where c and G are the values of the ve¬ 
locity of light in vacuum and Newton’s gravitation con¬ 
stant. Moreover, the gravitational wave averaged over 


several cycles leads to (cos (cut)), which is equal to 1/2. 

- Third, the energy Eabs that is absorbed by the star per 
unit of time by each resonant mode N averaged for a few 
cycles, is equal to the product of the gravitational wave 
force Egw(t) and the velocity of the mode £,m where £,m is 
the amplitude of the stimulated quadrupole mode. The 
gravitational force Egw(t) is equal to A/LMLnlim, where 
A/Lm and are the modal mass and the modal length of 
the mode, both of which depend on the properties of the 
acoustic eigenfunction, and 11^ is the second derivative 
of m’s spherical component, the perturbation tensor Hq. 
The modal length of the mode Ln is an effective length of 
the mode equivalent to the size of a resonant detector. Ln 
is unique for each acoustic mode, = 1/2 R Ixnl where 
Xn is a coefficient that depends of the eigenfunction £,m : 


Xn = 


dTtp^ 


p(t) [£,r,n(r) + 3£,h,n(i')] Udt. (4) 


where p and p^ are the density profile and mean density 
of the star, and E^r,n and £,H,rL are the radial and horizon¬ 
tal components of the quadrupole modal eigenfunction 
£,m , respectively. As originally computed for a resonant 
detector (e.g., Maggiore 2008), the absorbed energy per 
unit time is expressed as 


dEcLbs 

dt 


(Fgw(t)tN) 


(cu^ — 0)2 ]2 _|_ 4-p|2 (Jj2 ’ 
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where r\n is the damping rate of the mode. 

Equations (3) and (5) relate to the energy of the in¬ 
cident gravitational wave with the energy absorbed by 
the acoustic mode N, from which is obtained the average 
absorption cross-section: 


^abs(^) 


MlMTIrLh^ CU^ 327 rG 
(cu^ — -h dri^cu^ 


(6) 


The response of a mode to a gravitational perturba¬ 
tion is better evaluated by the integral of the absorp¬ 
tion cross-section, = J crabs(^^)da)/27r for which the 
limits of the integral correspond to the frequency inter¬ 
val of the gravitational wave-packet. However, as the 
function crabs(^^) only is significantly different from zero 
near each resonance frequency 0)^, conveniently, without 
much loss of accuracy, the limits of the integral can be 
replaced by —oo and -foo. For the same reason, the sec¬ 
ond term in the right-side of the equation (6) is expanded 
in (rin/cun). 

Ln near the resonance frequency, reduces to 


=7tXn 


M.]m G 
c 



(7) 


thus, Misi = Him M where Em is the normalized iner¬ 
tia of the mode. This approximate result shows that 
the integrated cross-section is independent of the quality 
factor of the mode, Qn = a)rL/(2r|rL)- This occurs be¬ 
cause at the peak cr(cu) is proportional to Q^l, so that 
Qrt cancels in J dcucr(cu). Using equation (1), Ln = 
Ert (/VLmG/c)(Rcu/c)^ where E^ = nXn (n-h 1 + ao)^. In 
the particular case that the sound speed is constant, for 
instance equal to the average sound speed Vg, (b = ttVs/R. 
Ln becomes identical to the one computed for a reso¬ 
nant mass detector, Ln = E* (MmG/c) (Vg/c)^ where 
E^ = Tt^Xn (n -I- 1 -f- ao)^ (Maggiore 2008). 

The relevant observable that will allow us to iden¬ 
tify the impact of gravitational waves in an oscillation 
mode is the phostopheric velocity that is computed as 
V^(cUrL) = l/(2r|Ti^dM) dEabs/dt (Goldreich & Keeley 
1977). Therefore, using equation (5), we obtain 


Vn(cU^) 


lT.*l-na)2 




(8) 


In the previous equation, we introduce the numerical con¬ 
stant oCs which is equal to 2 a/ 2, for which the exact value 
is fixed by observations. Since the amplitude of the strain 
decreases with the inverse of the distance, is computed 
as = de / d^ , where and d* are the distances of 
the source of gravitational radiation to the Earth and to 
the Sun-like star, and is the current strain prediction 
for the Earth detectors. An illustrative example is shown 
in Figure 1. 

The photospheric velocity, as shown in equation (8), 
relates to the excitation of quadrupole modes of fre¬ 
quency cUrt by a monochromatic gravitational wave with 
the same frequency cu, i.e., cu = 0)^. Nevertheless, in 
some cases as noticed by McKernan et al. (2014), cu drifts 
across the frequency uOn like it occurs during the inspiral 
phase of binary systems. The magnitude of the frequency 
drift cu is mostly dependent of the chip mass of the bi¬ 
nary Me, as cu 3 (GMc/c^)^^^ cu^^/^ (Maggiore 2008). 


The contribution of cu for the excitation of a quadrupole 
mode is only relevant if the duration of the gravitational 
time Tf = 1/v^ is smaller than the damping time of the 
mode Trt = l/[2r[n), or if the ratio T* = Tf/Xn is smaller 
than one. 

The steady-state (or saturated) limit corresponds to 
T^ > 1 (or Tf > Trt) for which cu can be neglected in 
the calculation of the photospheric velocity (equation 8). 
However, as shown by McKernan et al. (2014) in the re¬ 
verse case of undamped oscillations, for which T* ^ 1 
(or Tf ^ Trt), Iho slow frequency variation of the grav¬ 
itational wave must be taken into account, accordingly 
cUrt must be replaced by LVn + cUrtT^, where t is the time 
difference to the resonance. The calculation of Tf is given 
by 4.81 X 10^ (Mc/Mq)~^^^ (cuTv/mHz)”^^^^ (Maggiore 
2008). 

Rathore et al. (2005) have estimated that the averaged 
energy transfer to quadrupole mode by a slow varying 
frequency gravitational source is given by n}\^/[dd)n)- 
Accordingly, the energy absorbed by a stellar quadrupole 
mode per unit time in the undamped limit (T* ^ 1, sub¬ 
script u) is expressed as (dEabs/dt)^ = T* (dEabs/dt) 
where (dEabs/dt) corresponds to the solution given by 
equation (5). (dEabs/dt)^ relates with the work done by 
an external slow varying gravitational source (see McK¬ 
ernan et al. (2014) for details). Accordingly, the pho¬ 
tospheric velocity in the undamped limit is reduced in 
relation to the steady-state limit (i.e., V^t), it reads 

Vn.ultWn] = Vn(cUn)- Equally, CTabs,u(<^) and 

LrL,u are also reduced in the case of the undamped limit, 
as these quantities are obtained by multiplying crabs(^^) 
(equation 6) and Ln (equation 7) by the same T* factor. 

The integrated absorption-cross section (equation 7) 
and the photospheric velocity (equation 8) for a star can 
be scaled to the Sun. Using the eigenfrequency equation 

(1) retaining only the first term, and the scale relation 

(2) , the integrated absorption-cross section reads 



and the photospheric velocity reads 



where the subscript © denotes the Sun. The equiv¬ 
alent expressions for the undamped/unsaturated case 
(U ^ 1) are obtained by multiplying the right-side 
of both previous equations by the ratio (T*/T®) and 
(T*/T®)^/^, respectively. 

4. THE SENSITIVITY OF STAR’S DETECTORS TO 
GRAVITATIONAL WAVES 

Similarly to convectional resonant mass detectors, the 
stimulation of quadrupole oscillations in stars by an ex¬ 
ternal gravitational radiation source depends on the inte¬ 
grated cross-section and the quality factor. Table 1 lists 
these quantities, as well as the amplitude of the photo¬ 
spheric velocity for the quadrupole acoustic modes in the 
Sun. Moreover, we also computed the photospheric ve¬ 
locity in the case of undamped oscillations and the ratio 
T*. In this study we chose as fiducial gravitational wave 
source the AM CVn binary system (Me = 0.248 M©), 
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Fig. 1.— Impact of a gravitational wave in a Sun-like star detector. The photospheric velocity VrL(cUrL) (in cm s for quadrupole 
acoustic modes is a function of the distance to the gravitational wave source, as predicted by equation (8). In the case of modes with n ^ 5 
(see Table I), the photospheric velocity is computed in the undamped limit, i.e., the photospheric velocity corresponds to Vn,u = Vn- 

As an example of a gravitational wave source, we choose the compact binary system AM CVn (e.g., Solheim et al. 1998; Roelofs et al. 2007) 
located at a distance d0 = 606 pc, which has a size smaller than 0.5 of the solar radius and a orbital period of ~ 1029s, for which the strain 
amplitude for Earth detectors is predicted to be Ii 0 = 2 10“^^. Accordingly, the strain amplitude H* of the arriving gravitational wave to 
the Sun-like star is computed as = 3 / d^[A\l) . The vertical dashed-lines correspond to the distance of the gravitational wave 

source to the star’s detector: Sun (magenta) and a Sun-like star located at a distance of I AU (red) or 10, 100 and 1000 AU (blue). The 
vertical green and yellow areas correspond to 0.5 and 1.0 of the solar radius. The oblique lines correspond to the quadrupole modes: pt 
with 1 = 0,1,2,4, 5,6 (red), pt with i = 3, 7,12,18 (blue, green, magenta, cyan) correspond to H* = 3 10“^^/ d*(AU) with H 0 = 10“^^ 
(see Table 1). For reference, we show po and pi (black lines) corresponding to the steady-state (or saturated) limit. The horizontal cyan 
band (Vn ^ 1.0 cm s~^) defines the expected Vn threshold of the Plato mission. To illustrate the variety of sun-like star oscillations, we 
show several maximum Vn measurements (horizontal lines): Sun (green), Procyon (blue) and v Ind (magenta) stars (see main text). See 
Table 1 for further information about quadrupole modes. 


for which Tf is equal to ~ 1.54 x 10^^ sec. 

Existing observations and other theoretical quantities 
computed from an updated standard solar model (e.g., 
Lopes & Turck-Chieze 2013) are also summarized in the 
table. The quadrupole modes that have the largest inte¬ 
grated cross-sections and quality factors, as well as the 
largest photospheric velocity, are the modes of lower or¬ 
der. In the following, we give the reasons why Sun-like 
stars are good gravitational wave detectors: 

- It follows from equation (7) that the integrated cross- 
section of a mode depends on the modal mass and the 
length of the mode, as well as the sound speed in the in¬ 
terior of the star. Because stars have much larger masses 
than convectional resonant mass detectors, their inte¬ 
grated cross-section is many orders of magnitude larger. 
In the Sun, low order modes have a Ln that is 17 or 
20 orders of magnitude (if we consider the steady-state 
or undamped limits respectively) larger than the cur¬ 
rent detectors, such as Mario Schenberg (Aguiar et al. 
2006), Minigrial (Gottardi 2007) and Allegro (Mauceli 


et al. 1996). Gottardi (2007) estimates a Ln of the or¬ 
der of 9.8 10“^^cm“^ Hz for a high performance resonant 
mass detector. Even for high order modes (n ~ 15) the 
Sun, as a detector, performs better by a ten-magnitude 
factor. More massive stars have a fractionally smaller Ln 
than the Sun. Hence, from equation (9), it follows that 
Ln increases with M.^ and decreases with R, as for most 
of the observed Sun-like stars, like red giants (Mosser 
et al. 2012), the smaller mass corresponds to 0.5 Mq and 
the largest radius to 15 R©, which leads to a maximum 
reduction of Ln by a negligible factor of 0.02. The vari¬ 
ation of Xtl and Ln among these stars also introduces 
important corrections to Ln, in many cases, this should 
increase the value of Ln- Therefore the net result found 
for the Sun also holds for many of these stars. 

- Stars respond to gravitational radiation as a truly 
high Q oscillator. As shown in table 1, this is particu¬ 
larly relevant in the case of low order modes. Actually, 
the Qrt for small n is 1 to 2 orders of magnitude bet¬ 
ter than for current detectors. Gottardi (2007) estimates 
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that Qrt 10^ for a CuAl alloy spherical detector. In 
Sun-like stars, such high Q is due to the predicted low 
rates, believed to be caused by the damping processes 
related with radiation and turbulence convection of the 
upper layers of the star. This behaviour is most relevant 
for main sequence and subgiant stars (with A/L varying 
from 0.95 A/Lq to 2.0 A/l©) for which current theoreti¬ 
cal models predict that the rirt foi* low n modes is of 
the order of 10“^ to 10“^ (Houdek et al. 1999; Baudin 
et al. 2005). It was also found that for global acoustic 
modes (1^3), that r\n is independent of I and always 
increases with n. Besides, rirt always decreases with the 
mass of the star and also as its age increases. This be¬ 
haviour is similar for red giant stars. Recent seismic 
data has shown that the average damping rate (tIti) of 
Sun-like stars increases with the effective temperature. 
This relation holds not only for main sequence and sub¬ 
giant field stars, but equally for red giant stars, some of 
which were discovered in open stellar clusters (e.g., Cor- 
saro et al. 2012). As the (rirt) of red giant stars is two 
orders of magnitude smaller than for the Sun’s case, we 
expect that the rirt of low order modes varies by a sim¬ 
ilar amount, leading to an increase of Q^. Nonetheless, 
a smaller quality factor could be beneficial, as a lower 
Q factor increases the possibility of quadrupole modes 
being easily excited by gravitational radiation. 

- According to equation (8), the low order quadrupole 
modes are the ones for which the photospheric velocity 
is predicted to have the largest value. The velocities of 
the modes of higher order are not so high. The increase 
of the photospheric velocity depends of the quantities 
on the ratio a)^/r\n which for low n have not yet been 
measured for the Sun, and only theoretical predictions 
are available, as shown in italics in table 1. Neverthe¬ 
less, both quantities are predicted by theoretical models, 
which successfully reproduce the observational data at 
higher frequencies. The Ln and its related quantity Xn 
are the other critical quantities that affect the photo- 
spheric velocity. Unlike for a typical detector, Ln and 
Xn decrease rapidly with n. As shown in equation (4), 
this is due to the fact that the density in the star de¬ 
creases with the increase of the radius, and the depth 
reached by the acoustic eigenfunction £,m decreases as n 
increases. As shown in table 1, Ln is only significant for 
modes of low order. The decrease of Ln and r|k^ with n 
is the main reason for the rapid decrease of photospheric 
velocity with n. Nonetheless, we notice that in the case 
of the frequency of the gravitational wave drifts during 
the excitation of the quadrupole mode, the maximum 
attained will be reduced by up to several orders of magni¬ 
tude, as we discuss in the previous section. Table 1 shows 
the photospheric velocity computed for a strain of 10“^^, 
which for the first low order modes, is predicted to have a 
maximum amplitude of 10“^ cms“^ or 7 10^^ cms“^ re¬ 
spectively for the steady-state (T* ^ 1) and undamped 
(T* l) limit cases. In this study, the reduction of 
over steady-state is as much as two orders of magnitude, 
and will be greatest for modes with the largest cross- 
sections. Indeed, the excitation of quadrupole modes by 
gravitational waves is limited at low order by the fact 
that oscillations are unsaturated, and at higher order 
modes because these have small cross sections to the im¬ 
pact of gravitational waves (see Table 1). 


Thus in order for gravitational waves in the star acous¬ 
tic spectrum to be detected, it would be necessary for 
to be larger than the current signal-yo-noise ratio instru¬ 
mental threshold. A current estimation of this threshold 
for a future interplanetary space mission (Appourchaux 
et al. 2009), made from an average of the 50 modes ob¬ 
served by the GOLF experiment (Turck-Chieze et al. 
2004; Garcia et al. 2007; Jimenez & Garcia 2009) dur¬ 
ing a period of 10 years fixes the observation limit of 1 cr 
level between 10“^cm s“^ and 3 10“^cm s“^. Neverthe¬ 
less, for the next planed seismic missions, TESS (Ricker 
et al. 2014) and Plato (Rauer 2013) the minimum Vn 
measured is expected to be of the order of 1 cms“^. 
Although, this precision is not sufficient to detect grav¬ 
itational waves with current sun-like stars, it is possible 
that goal to be achieved in a next generation of seis¬ 
mic instruments, as we show in figure 1 if a precision of 
10“^ cm s“^ is attained, this type of detection could be 
successful in some specific multi-body stellar systems. In 
the remaining of the article, this value is used as a de¬ 
tectability threshold for a future space mission. Indeed, 
it follows that the Vn for the f-mode, pi-mode and p 2 - 
mode increases above the threshold of 10“^ cms“^, if 
d* is smaller than 10 AU, as illustrated schematically in 
Figure 1 for the case of a fiducial compact binary as the 
AM GVn binary. Obviously, if the strain is 2 orders of 
magnitude larger, it will be much easier to successfully 
detect the gravitational radiation (cf. Figure 1). 

The variation of mass and radius of the star affects 
the amplitude of the photospheric velocity, and it fol¬ 
lows from equation (10) that Vn increases with A4 and 
decreases with . Thus ignoring the variation of Xn and 
Tin, and taking the minimum A4 and maximum R val¬ 
ues from the recent list of observed Sun-like stars (e.g., 
Mosser et al. 2012; Ghaplin et al. 2014), at 0.5 Ad© and 
15 Rq, the photospheric velocity is reduced by a factor 
2 10“^. One such example is the star KIC 5822889 for 
which Vn is O.OlVn© For other stars the effect is the re¬ 
verse, as for KIC 7970740 for which Vn ~ 1.3Vn0 (Chap¬ 
lin et al. 2014). This is a rough estimate, because these 
calculations do not account for the potential differences 
coming from T*, r|n and Xn- la particular, notice the 
case of quadrupole modes with very low n for which the 
excitation by gravitational waves occurs in the unsatu¬ 
rated case, T* = Tf/Tn reduces further the value of Vn- 
In the case of the AM GVn compact binary, T* ~ 10“^ 
assuming that r|n is similar to the Sun, Vn reduces by 
a further factor. For instance, r|n for low n is of 

the order of 10“^pHz and does not change much among 
this type of star, as Xn increases for more massive stars, 
due to the increase of the density in the star’s core. Al¬ 
though Sun-like stars can have different Xn and rjn, the 
low n modes should be always the ones with the largest 
velocity amplitudes. 

Actually, in many spectra of Sun-like stars, acous¬ 
tic modes were found with large photospheric veloci¬ 
ties, believed to be stimulated by the turbulent convec¬ 
tion in the envelopes of these stars. Figure 1 shows the 
maximum photospheric velocities measured for a few of 
the most well-known stars ^ other than the Sun - with 

Characteristics of a few stars: Procyon has M = 1.5M.Q and 
R = 2.OR©; and v Ind has M = O.SSM.© and R = 2.OR© 
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20 cms ^ (Jimenez & Garcia 2009), Procyon with 38 
cms“^ (Kervella et al. 2004; Leccia et al. 2007; Arentoft 
et al. 2008; Bedding et al. 2010) and the sub-giant star y 
Ind with 650 cm (Bedding et al. 1996, 2006; Carrier 
et al. 2007; Kjeldsen et al. 2008). This behavior fol¬ 
lows the well-known scaling equation for the maximum 
velocity amplitude of acoustic modes Vn,max excited by 
turbulent convection, which increases with the luminos¬ 
ity and decreases with mass and effective temperature of 
the star (Kjeldsen & Bedding 2011). This equation re¬ 
mains valid for red giant stars, as found in observations 
for which Vn,max varies from 10 to 800 cm s~^ (Samadi 
et al. 2013). ’ 

The acoustic mode frequencies of Sun-like stars can 
probe gravitational waves in the range from 10“^ Hz to 
10“^ Hz, which overlaps with the gravitational radiation 
frequency range that will be probed by eLISA (Amaro- 
Seoane P. et al. 2013) for the high frequency range and 
EPTA (Ferdman et al. 2010) and SKA (Johnston et al. 
2007) in the low frequency range. More significant even is 
the possibility of probing the range from 10“^ Hz to 10“^ 
Hz, a frequency range for which no experiment has yet 
been planned. This is illustrated in Figure 2. In the low 
frequency range, this region corresponds to the predic¬ 
tions of stochastic background radiation and supermas- 
sive binaries. In the high frequency range, it corresponds 
to unresolvable galactic binaries, extreme mass ratio in¬ 
spirals and resolved galactic binaries (Sathyaprakash & 
Schutz 2009). One of the targets of eLISA will be nearby 
ultracompact binaries, such as the binary system AM 
CVn discussed in this article as a template. The strat¬ 
egy of using Sun-like stars as detectors will enable us to 
determine the impact of gravitational waves in the pho- 
tospheric velocities of quadrupole modes, or, at least, to 
fix an upper limit on the strain of the gravitational waves 
generated by these sources. The h sensitivity of eLISA 
will be 10“^^ to 10“^^ for the frequency range of 10“^ 
Hz to 10“^ Hz which is not sufficient to detect the gravi¬ 
tational waves produced by the AM CVn binary system, 
predicted to have a strain amplitude on Earth of 10“^^ 
(cf. Figure 2). Yet, a future seismic mission could detect 
such gravitational waves in stars like the Sun located at 
a distance of either 10 AU or 1000 AU from this binary 
system, if the is of the order of 10“^^ or 10“^^, re¬ 
spectively. As shown in figure 1 the quadrupole modes 
of lowest orders will be stimulated by these gravitational 
waves producing photospheric velocities possibly above 
the observational limit of the detector. Moreover, other 
stars like sub-giant and red giant stars could scan other 
parts of the frequency range of the gravitational wave 
spectrum, including outside of the current range of de¬ 
tectors. Red giant stars with oscillations within the fre¬ 
quency range from 10“^ Hz to 10“^ Hz (Mosser et al. 
2013) can be used to explore events related with super- 
massive binaries for which the strain on Earth is pre¬ 
dicted to be of the order of 10“^^ (see Moore 2014, and 
references therein). 

5. STELLAR SYSTEMS AND THE DETECTABILITY OF 
GRAVITATIONAL WAVES 

Central to the applicability of using sun-like stars as 
detectors is to find an oscillating star for which the pho¬ 
tospheric velocity is above the threshold of detectability 
(cf. Figure 1). As discussed previously for the sun’s case. 


the photospheric velocity is of the order 10“^ cm s“^ for 
a strain of h_ 2 o = 1 (see Table 1). Therefore for such 
a signal to be detected it is necessary to find a way to 
increase the photospheric velocity by 4 orders of magni¬ 
tude. Two obvious possibilities are to choose stars for 
which their Xn sensitive to gravitational waves is larger 
than the sun (like red-giants), or stronger gravitational 
wave sources (like the coalescence of black holes with a 
typical strain ~ 10~^^). Here, nevertheless the discussion 
is focused in another point: the possibility of the next 
generation of stellar missions discovering gravitationally 
bound multi-body systems with a compact binary and 
an oscillating sun-like star, for which the seismic instru¬ 
ment has the necessary accuracy to measure the impact 
of gravitational waves in the acoustic oscillations. 

The probability of finding an oscillating star nearly 
enough to a compact binary is quite small, nevertheless, 
recently many gravitationally bound multi-body systems 
of three or more stars have been discovered, which gives 
some hope that one of such unique systems could actually 
be found. The current theory of stellar formation, pre¬ 
dicts the existence of many binary, triple and higher mul¬ 
tiplicity stellar systems. This has also been confirmed by 
many astronomical observations. Tokovinin (2014) stud¬ 
ied the multiplicity of multi-body systems within 67 pc 
of the Sun and found that 13% of stellar systems have 
three or more components. Raghavan et al. (2010) esti¬ 
mated that as much as 8% are three-body systems, and 
Rappaport et al. (2013) using Kepler data estimated that 
20% of close binaries have tertiary companions. More¬ 
over, Holberg et al. (2002) have reported that within the 
first 20pc of the Sun, 25% of white dwarf are in binary 
systems. 

Several triple systems have been discovered with char¬ 
acteristics relatively close to the ideal case discussed in 
this article, compact binaries (Meliani et al. 2000) with a 
star companion. In particular, the COROT and Kepler 
missions have discovered several specific triple star sys¬ 
tems (Southworth 2014), like the KOI-126, a very small 
hierarchical triple system composed by a short-period bi¬ 
nary of two 0.2M© stars orbiting by an evolved G-star 
(with a mass of l.fiM© and radius of 2.0 R©) (Carter 
et al. 2011). The semi-major axis of the inner and 
outer binaries are 0.021 AU and 0.24 AU respectively. 
Another equally different example is the triple system 
J0337+1715, a compact system with a total dimension 
smaller than 2 AU, where a close binary of a white dwarf 
and neutron star is orbited by a second white dwarf (Ran¬ 
som et al. 2014). Equally, Kilic et al. (2014) have discov¬ 
ered a close binary of white dwarfs WD 0931-1-444 that 
possibly have an M dwarf companion, although for the 
moment sun-like oscillations have not yet been discovered 
in M dwarfs (Rodrfguez-Lopez et al. 2015). 

Finally there is the possibility, that among the many 
binary systems formed by a white dwarf and sun-like 
star discovered by the Kepler mission (Rappaport et al. 
2015; Southworth 2014), like KOI-3278 (Kruse & Agol 
2014), the white dwarf could actually be a close compact 
binary. Equally, higher order multiplicity system could 
also be worth investigating to look for compact binaries 
and oscillating stars (Riddle et al. 2015), like some 2-|-2 
quadruple systems which have an outer separation of 500 
AU. 


Sun-like star's detectors 



Fig. 2 .— Comparison of Sun-like star’s detectors with current gravitational waves detectors: (i) The range of frequencies for probing 
for gravitational wave sources is evaluated from the observed acoustic spectra of Sun-like stars: main sequence and subgiant stars show 
oscillations from 2 10“"^ Hz to 10“^ Hz (Arentoft et al. 2008); and red giant stars from 10“^ Hz to 5 10“^ Hz (Mosser et al. 2013) (magenta 
shadowed area), (ii) The range of strain amplitudes possible to probe by natural star detectors^ depends mainly from the amplitude 
generated by the gravitational wave source, as unlike for gravitational wave detectors on Earth, Sun-like stars can be located quite close 
to the source. As for the example shown in Figure 1 if a natural star detector is located at a distance of 1 AU from the gravitational wave 
source (red vertical line in Figure 1), for which H*(1AU) ~ 3 10“^"^, this corresponds to H 0 ~ 210“^^ on Earth, leading to an effective 
strain gain of 10+®. Equally, if the star detector is located at 1 solar radius from the gravitational wave source H^( 1 R 0 ) ~ 510“^^. At 
these distances all the low order modes (n ^ 4) of the star are stimulated by the gravitational waves generated by the compact stellar 
binary. The strain sensitivity curve of the different gravitational detectors and most common sources of gravitational waves were adapted 
from Moore (2014). 


Although asteroseismology is not as precise as helio- 
seismology, if the optimal triple or higher order multi¬ 
body stellar system is found the accuracy of the future 
(and possibly present) stellar seismic missions could be 
sufficient to detect gravitational waves. Currently, the 
Plato mission observations for quiet stars, expects that 
for an exposure time of only 15 minutes to be sufficient 
to average out the perturbing signal well below 1 m s“^, 
and in some cases the noise even gets down 10 cm s“^ in 
about 20-30 minutes (Rauer 2013). It is possible that in 
cases of much longer time exposures, the compute pho- 
tospheric velocities will be well below 1 cm s~^ (Rauer 
et al. 2014). 

Altogether, by using this proposed method, main- 
sequence, sub-giant and red giant stars increase signif¬ 
icantly the detection power and the spectral range for 
identifying gravitational wave imprints in acoustic oscil¬ 
lations. 

6 . SUMMARY AND CONCLUSION 

Stars like the Sun have the potential to become ideal 
detectors of gravitational waves. Not only do they have a 
high scattering cross-section for gravitational waves, high 
quality factors and very likely large photospheric veloc¬ 
ities, but they also may be located near strong gravita¬ 
tional wave candidates. Helio and asteroseismology allow 


us to probe these natural detectors with very high preci¬ 
sion. Indeed, for each oscillating star, the radial modes 
and dipole modes can be used to determine the oscillator 
characteristics like the damping rates, which is critical to 
allow us to disentangle the stimulation caused by gravita¬ 
tional radiation from the excitation produced by turbu¬ 
lent convection. The standard theory of stellar pulsations 
by turbulent convection establishes that low order acous¬ 
tic modes (1 ^ 3) are equally excited and their excitation 
depends uniquely on the frequency of the mode and not 
on the degree. This model has already been proven to 
be consistent with observations for a large range of fre¬ 
quencies. Therefore, the radial and dipole modes with 
frequencies near the frequencies of the quadruple modes 
of low order can be used as reference to isolate the stim¬ 
ulation of low order quadrupole modes by gravitational 
waves from the stellar self-excitation. 

Although main-sequence, subgiant and red giant stars, 
all have the possibility to be good detectors of gravita¬ 
tional waves, subgiants and red giants have a substantial 
advantage in relation to main-sequence stars, as their rich 
spectrum has shown the existence of quadrupole mixed 
modes (Benomar et al. 2013). These modes behave like 
gravity waves in the star’s core and like acoustic waves 
in the stellar surface. Accordingly, their eigenfunctions 
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have large amplitude in both the core and the surface. 
As a consequence the modal length, is much larger 
for these stars than in the Sun’s case, so that their pho- 
tospheric velocities are larger than the ones predicted for 
the Sun. 

The Sun and similar stars also have quadrupole gravity 
modes that can be excited by gravitational waves. These 
modes have very large modal lengths Ln compared with 
acoustic modes, because eigenfunctions of gravity modes 
have larger amplitudes in the star’s core. Unfortunately, 
gravity modes are strongly attenuated when they prop¬ 
agate through the convection region and as such it is 
very difficult to observe them in the star’s surface. Actu¬ 
ally, in the Sun it has been a matter of dispute whether 
a few gravity modes candidates were observed success¬ 
fully (Turck-Chieze et al. 2004). 

The group of all oscillating stars in the solar neighbor¬ 
hood within one thousand parsecs radius, constitutes the 
largest detector ever for gravitational radiation. Stars 
have some advantages over the current Earth detec¬ 
tors. There are thousands of oscillating stars scattered 
throughout space, some of which can be found relatively 
near gravitational wave sources. Alignments of stars 
between the source and the Solar System can monitor 
the progression of gravitational waves throughout space, 
which can be use as a test to probe General Relativity, 
a goal that is difficult to achieve with present man-made 
detectors. 

Equally important is the fact that such a new method 
to prove gravitational wave radiation could probe an im¬ 
portant part of the spectral gravitational radiation win¬ 
dow - from one micro Hz to 100 micro Hz (cf. Eigure 2), 
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